Dominant gamma-ray bursts production in the early universe 
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ABSTRACT 

It has been known that at least some of the observed gamma-ray bursts 
(GRBs) are produced at cosmological distances and the GRB production rate 
may follow the star formation rate. We model the BATSE detected intensity 
distribution of long GRBs in order to determine their space density distribution 
and opening angle distribution. Our main results are: the lower and upper 
distance limits to the GRB production are z ~ 0.24 and z > 10, respectively; the 
GRB opening angle follows exponential distribution and the mean opening angle 
is about 0.03 radian; the peak luminosity appears a better standard candle than 
the total energy of a GRB. 

Subject headings: gamma rays: bursts - stars: formation - cosmology: early 
universe, reionization 



1. Introduction 



Recent research results indicate that radiation from gamma-ray bursts (GRBs) may be 
constrained mostly in narrow beams (Waxman et al. 1998; Fruchter et al. 1999; Woosley 2001), 
which suggests that GRBs may be standard candles at cosmological distances (Frail et al. 2001). 
Previous studies of the GRB space density distribution suggested that GRB production rate 
may follow the observed star formation rate for small values of redshift (Wijers et al. 1998), 
but may increase monotonically for very high values of redshift (Schaefer et al. 2001; Norris 2002) 
however neither the lower limit nor the upper limit of redshift for GRB production could be 
determined reliably previously because these studies are either limited by statistics or have to 
rely on some empirical relationships between some statistical properties of selected samples of 
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GRBs. Because both the space density distribution and beaming effect play important roles 
for the observed GRB intensity distribution, here we make use of the BATSE/CGRO GRB 
intensity distribution, which is the largest sample of GRBs collected so far, to determine 
directly the GRB space density and opening angle distribution. 



2. Models for GRB intensity distribution 

For the GRB space density distribution, we assume the following model, under the 
cosmological parameters of Q = 1, Hq = 50 km s _1 Mpc -1 and go — 0.5, in order to compare 
with previous star formation rate measurements (Steidel et al. 1999). 

r o, n(z < zo) 

n(l + z) = < n (l + zf- 5 , if(z Q < z < z break ) . (1) 

Zbreak ) 

The functional form is taken from the measured star formation rate (Steidel et al. 1999) 
with z > 0.2 and Zb re ak ~ 1- Our main goal in this paper is to determine Zq and Zbreak f° r 
GRBs from the observed GRB intensity distribution, in order to answer the question where 
in the universe GRBs are produced predominantly. 

For the GRB opening angle distribution, we study three cases when the GRB opening 
angle 9 follows Gaussian, exponential, or power-law distributions. For a Gaussian distribu- 
tion of 9, the probability that the opening angle equals to a given angle is, 

1 (e-»n) 2 

P{9) = e-^K (2) 
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We assume that the total energy E is the same for for every long GRB, i.e., they can be 
used as standard candles(Frail et al. 2001; Balazs et al. 2003) and the opening angle 9 is 
very small, the probability that the observed fluence / is higher than a specified value F is 
then, 
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where is the luminosity distance. Thus the number of bursts detected is, 



M, (3) 



N — -Trnd 2 L dd L / ^d9. (4) 
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Similarly, for an exponential distribution with the parameter A, 



N 



-nnd 2 L dd L 




6 2 \e- xe d9, 



(5) 



and for a power-law distribution with the parameter t, 



N = 




(6) 



If we assume the peak luminosity of GRBs is the standard candle, we only need to 
replace E and F by the standard candle luminosity L and the peak flux P in equations 
(3), (4), (5) and (6), respectively. We number these slightly modified equations as (3)', (4)', 
(5)' and (6)', respectively; for brevity we do not list these equations explicitly. 



We apply equation (1) for the GRB space density distribution and equations (4), (4)', 
(5), (5)', (6) and (6)' for the GRB opening angle distribution to fit the observed GRB intensity 
distribution. We focus our study on long GRBs (T 90 > 2s) (Kouveliotou et al. 1993) because 
currently only long bursts are known to be originated at cosmological distances from the 
direct redshift measurements of the optical afterglows of some long GRBs. All long GRBs 
with peak flux values and detection probability above 5% included in the GRB 4B+ Catalog 
(Paciesas et al. 2000) are selected, resulting in 1727 GRBs as our sample. The fluence values 
of these GRBs are taken also directly from the catalog. For the peak flux distribution, the 
detection efficiency correction is done directly using the detection probability of each peak 
flux value given in the catalog. Because the fluence and peak flux does not have a good linear 
correlation, for each fluence range we average the detection probability of all GRBs included 
in this fluence range and then correct for the detection efficiency of this fluence range. 



3. 



Results 
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Total energy as standard candle. 



Peak luminosity as standard candle. 





Gaussian exponential power-law 


"K.T 


281.3/194 261.1/195 623.2/195 




# =0.028±0.007 

n n^n^n nnc A=38.1±3.8 t=4.51±0.09 
(j=0.050±0.008 


z 

Zbreak 


0.23 ±0.07 0.23±0.08 0.27 ±0.11 
7.7 - oo 7.9-oo 7.2 - oo 


n 


231 ±51 215 ±45 207 ± 44 
0.69 ±0.37 0.80 ±0.41 0.72 ±0.36 




Gaussian exponential power-law 


'X.r 


241.1/194 227.1/195 539.7/195 




fl =0.031±0.008 

( x=0.047±0.01 ^=34-4±2.1 t=4.52±0.09 


z 

Zbreak 


0.23 ±0.09 0.24±0.07 0.29 ± 0.09 
9.5 - oo 9.8-oo 9.2 - oo 


n 
L 


273 ± 44 264 ± 49 257 ± 39 
0.67 ±0.32 0.64 ±0.22 0.69 ± 0.30 



Table 1: Xr = X 2 /dof. L is in units of 10 51 ergs/s. n in units of Gpc 3 /year. All errors 
quoted are for 68.3% confidence. For both cases the exponential distribution of the opening 
angle offers the best fits; however the residuals for the combination of peak luminosity as 
the standard candle and the exponential distribution is the only statistically acceptable one. 




Fig. 1. — Comparison of different models with data. Top panel: long GRBs with fluence as 
standard candle; Lower panel: long GRBs with peak luminosity as standard candle. 
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Our fitting results are shown in Table. 1 and fig.l. Because our fittings are done with the 
integral distributions, the numbers of GRBs in different bins are not completely independent. 
To address this problem, we carried out extensive bootstrap simulations. Taking the 1727 
GRBs as the seed distribution, we randomly generated a large numbers of peak flux or fluence 
distributions and fitted the sampled distributions in the same way. The best values and errors 
of all parameters are then calculated from these fittings; the bootstrapped parameter values 
and errors are consistent with that listed in Table 1. We therefore believe the correction 
between the different points does not have any significant impact to our results. From the 
fitting results, we conclude: 

• For long GRBs, despite the different fitting residuals of the two standard candle mod- 
els and three opening angle distributions assumed, the cosmological parameters of 
GRBs inferred are strikingly similar, indicating the robustness of the GRB cosmologi- 
cal model. 

• The power-law distribution for the GRB opening angles cannot fit the data and is thus 
rejected confidently. 

• The Gaussian distributions for the GRB opening angles fit the data only marginally, 
and may be rejected too. 

• For both standard candle models the exponential distribution for the GRB opening 
angles provides the best fit to the data. 

• The combination of the peak luminosity as standard candle and the exponential dis- 
tribution for the GRB opening angles is the only one with statistically acceptable fit 
to the data; we therefore accept this as the accepted model for long GRBs. Under this 
model the peak luminosity of a GRB is (0.64 ± 0.22) x 10 51 ergs/s, the exponential 
constant of the opening angle distribution is 34 ± 2 (radian -1 ), the lower bound to the 
GRB redshift is z = 0.24 ±0.07, and the upper break redshift is Zbreak > 9.8 (no lower 
limit to the break is detected because not enough sufficiently faint GRBs have been 
observed with BATSE). 

4. Discussion 

First we compare our study with previous studies of GRB cosmological models with 
BATSE intensity distribution (Fenimore et al. 1993; Fenimore and Bloom 1995) 
(Wickramasinghe et al. 1993; Wijers et al. 1998). All these studies concluded that the dimmest 
GRBs are at redshift around unity, probaly following the measured star formation rate which 
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peaks also at redshift around unity (Steidel et al. 1999). The main differences between our 
approach and previous studies are the following: 

• We take the GRB opening angle distribution as part of the cosmological model. This 
is very important, since the observed GRB flux depends strongly on the opening angle, 
as well as the redshift. 

• We separate the long and short GRBs as two different classes in our study. Our 
conclusions are only applicable to long GRBs. For short GRBs, no conclusion can be 
drawn due to the limited statistical quality of data. 

• One simplification of our study is that we did not make the energy spectral correction 
to the fluence or peak flux as a function of z (K-correction); we just take the fluence 
or peak flux values as listed in the catalog. Technically this correction is not straight- 
forward, because we can no longer uniquely convert the fluence or peak flux directly 
to redshift, because we believe that GRB radiation is beamed, and the opening angle 
follows a certain distribution. Therefore the observed fluence or peak flux is due to 
both the redshift and the opening angle of the GRB. The good fit of the assumed 
simple cosmological model with only several free-parameters indicates that this cor- 
rection may not impact our conclusions significantly. In fact, since the energy spectra 
of GRBs in the BATSE band follows approximately a power- law with an index of —2, 
the K-correction is not important. 

Frail et al. (2001) suggested that the total energy of GRBs may be a constant, based 
on the opening angles of several GRBs determined from their afterglow light curve break. 
Our result is not in major conflict with their suggestion, but favors the peak luminosity 
as the standard candle. Perhaps the direction of the GRB beam is not completely stable, 
as indicated from observations of other relativistic jets and outflows from stellar mass and 
supermassive black holes; only when the beam is pointed directly at us can we observe 
the maximum flux from the GRB. However it has been argued that for some GRBs the 
lower limit of the total energy released is 10 52 ergs as determined from the redshifted X-ray 
emission lines (Ghisellini et al. 2002). However the reality of the detection of the emission 
lines is still under debate, and an alternative interpretation of the detected emission lines 
also exists (Wang et al. 2002). 



- 8 - 



10| - l l l l l l l 




Opening angle (radian) 

Fig. 2. — The data for the obsered GRB opening angles are taken from Bloom et al. 2003; 
arrows indicate lower and upper limits of GRB opening angles. The power-law distribution 
fits the data for large opening angles, but fails to predict the absence of GRBs detected 
with opening angles smaller than 0.04 radian. The Gaussian distribution agrees with the 
peak reasonably well, but cannot predict enough GRBs with opening angles greater than 0.2 
radian. The exponential distribution is consistent with the data. 
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The opening angles of some GRBs with afterglow detections have been determined 
previously from the GRB afterglow data (Frail et al. 2001; Bloom et al. 2003); we call this 
Frail-distribution. In fig. 2 we compare the Frail-distribution with the opening angle distribu- 
tions determined with the observed GRB intensity distribution data, just as a sanity check. 
The power-law distribution agrees with the Frail-distribution for large opening angles, but 
fails to predict the absence of GRBs detected with opening angles smaller than 0.04 radian. 
No GRBs with opening angles narrower than 0.04 radian was directly observed, but the 
recent report on the polarization of GRB021206 suggests that it is possible for some GRBs 
to have opening angles narrower than 0.01 radian (Coburn and Boggs 2003; Waxman 2003). 
The Gaussian distribution agrees with the peak of the Frail-distribution distribution reason- 
ably well; however this model cannot predict enough GRBs with opening angles greater than 
0.2 radian. The exponential distribution is consistent with the Frail-distribution distribution 
for both standard candle models, but favors slightly the peak luminosity standard candle 
model. It was suggested that for the model in which the GRBs follow the star formation rate 
strictly, the most possible opening angle should be around 0.1 radian (Bloom et al. 2003). 
However, if the majority of faint GRBs originate at higher redshift than the peak of star 
formation rate as our model predicts, then a narrower average opening angle is required. 

We therefore conclude that GRB opening angles follow an exponential distribution with 
a mean opening angle of 1/A = 0.03 radian, and the GRB space density follows a power- law 
with an index of —3.5 between z = 0.24 ± 0.07 and z break > 10. The general trend of 
the GRB space density distribution determined here as a function of z is consistent with 
previous studies (Schaefer et al. 2001; Norris 2002); however our result is statistically much 
more robust, and in particular zq is determined for the first time and the lower limit to Zbreak 
is constrained more strictly than before. Our result therefore demonstrates that GRBs are 
no longer produced frequently in the nearby and present universe at z < 0.24 ± 0.07, and 
GRBs are predominantly produced in the early universe at z > 10. Indeed most GRBs 
with direct redshift measurements have redshifts z > 0.3 (Djorgovski et al. 2001), with the 
exceptions of GRB 980425 (,2=0.0085) and GRB 030329 (z=0.1687) which may be associated 
with two peculiar but similar supernovae SN1998bw (van Paradijs 1999), and SN2003dh 
(Stanek et al. 2003). It remains to be seen if these GRBs belong to a subset of GRBs which 
are different from the majority of GRBs which seem to be produced at much larger redshift. 

It has been suggested that GRBs follow star formation rate, because it is commonly be- 
lieved that GRBs are produced as the final gravitational collapse of massive stars (Djorgovski et al. 2001; 
Paczynski 1998). The GRB space density determined here follows the measured star forma- 
tion rate between z ~ 0.24 to z ~ 1. The determination of z = 0.24 ± 0.07 suggests that 
the massive stars capable of producing the majority of long GRBs are no longer formed 
frequently today. Future GRB experiments capable of measuring the redshift of many more 
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GRBs and thus providing direct distance measurements of GRBs, such as the Swift mission 
to be launched in late 2003 (Gehrels 2000), will test this prediction firmly. 

The monotonic increase of the GRB space density up to at least z — 10 suggests 
that the formation rate of massive stars capable of producing GRBs should peak at sig- 
nificantly higher z than the measured star formation rate. This implies that the "dark 
ages" of the universe ended much earlier than previously believed at z ~ 4 — 5 based on 
the measured space density distribution of high-redshift QSOs (Hook et al. 2002). There 
are also some previous independent works showing that an essential amount of long GRBs 
should be at z > 4 (Choudhury and Srianand 2002; Bagoly et al. 2003) or even up to 
z = 20(Meszaros and Meszaros 1996). A recently detected GRB031203 by INTEGRAL 
was reported to be located at z > 9 (formally z > 11) (GRBlog 2003), consistent with our 
model prediction. Therefore GRBs may be a probe of the epoch of reionization of the early 
universe (Lamb and Haiman 2003). Currently the GRB intensity distribution data cannot 
constrain the exact upper limit of z break , because not enough extremely faint GRBs have 
been observed. Future GRB experiments, such as the proposed "Next Generation GRB 
Mission" (Grindlay et al. 2001) would detect many more GRBs originated in the extremely 
early universe, which may determine the exact ending time of the "dark ages" of the universe 
in order to test the prediction of some numerical simulations (Abel et al. 2002) which sug- 
gests that the first stars may have been formed at z ~ 20. GRBs may prove to be the most 
sensitive and ultimate probes into the early epoch of the universe when the first structures 
were formed. 

Acknowledgement: We thank Drs. Zigao Dai, Tan Lu, Yongfeng Huang and Darning Wei 
for valuable comments to the manuscript. This study is supported in part by the Special 
Funds for Major State Basic Research Projects (10233010) and by the National Natural 
Science Foundation of China. SNZ also acknowledges supports by NASA's Marshall Space 
Flight Center and through NASA's Long Term Space Astrophysics Program. 

REFERENCES 

Abel, T., Bryan, G. L., k Norman, M. L. 2002 Science, 295, 93 

Bagoly, Z., Csabai, A., Meszaros, A., Meszaros, P., Horvath, L, Balazs, L. G., & Vavrek, R. 
2003, A&A, 398, 919 

Balazs, L. G., Bagoly, Z., Horvath, L, Meszaros, A. & Meszaros, P. 2003, A&A, 401, 129 
Bloom, J. S., Frail,D. A. & Kulkarni, S. R. 2003, ApJ, 594, 674 



- 11 - 



Choudhury, T. R, & Srianand, R. 2002, MNRAS, 336, L27 
Coburn, w., & Boggs, S. E. 2003, Nature, 423, 415 

Djorgovski, S., et al. 2001, ESO Astrophysics Symposia, Berlin: Springer Verlag 

Fenimore, E.E., & Bloom, J.S. 1995, ApJ, 453, 25 

Fenimore, E.E., et al. 1993, Nature, 366, 40 

Frail, D. A. 2001, et al, ApJ, 562, L55 

Fruchter, A. S., et al. 1999, ApJ, 519, L13 

Gehrels, N. A. 2000, Procspie 4140, 42 

Ghisellini, G., et al. 2002, A&A, 389, L33 

http://grad40.as.utexas.edu/grblog.php?GCN=2468 

Grindlay, J., et al. 2001, Presented at Woods Hole GRB Conf., to appear in AIP Conf. Proc. 

Hook, I. M., et al. 2002, A&A, 391, 509 

Kouveliotou, C, et al. 1993, ApJ, 413, L101 

Lamb, D.Q., & Haiman, Z. 2003, AAS, HEAD2003, 16.15 

Meegan, C. A., et al. 1992, Nature 355, 143 

Meszaros, A., & Meszaros, P. 1996, ApJ, 466, 29 

Norris, J. P. 2002, ApJ, 579, 386 

Paciesas, W. S., et al, VizieR On-line Data Catalog: IX/20A. Originally published in: 2000, 
ApJ, 122, 465P; ApJ, 122, 497H 

Paczynsky, B. 1998, ApJ, 494, L45 

Schaefer, B. E., Deng, M. & Band, D.L. 2001, ApJ, 563, L123 
Stanek, K.Z., et al. 2003, ApJ, 591, L17 
Steidel, C. C, et al. 1999, ApJ, 519, 1 
van Paradijs, J. 1999, Science, 286, 693 



- 12 - 

Wang, W., Zhao, Y., & You, J. H. 2002, ApJ, 576, L37 
Waxman, E., Kulkarni, S. R. & Frail, D. A. 1998, ApJ, 497, 288 
Waxman, E. 2003, Nature, 423, 388 
Wijers, R.A.M.J., et al. 1998, MNRAS, 294, L13 
Wickramasinghe, W.A.D.T., et al. 1993, ApJ, 411, L55 
Woosley, S. 2001, Nature, 414, 853 

Correspondence should be addressed to S.N.Z. (e-mail: zhangsn@tsinghua.edu.cn). 



This preprint was prepared with the AAS IATf^X macros v5.2. 



